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We study the long time asymptotics of probability density functions (pdfs) of Levy flights in dif- 
ferent confining potentials. For that we use two models: Langevin - driven and (Levy - Schrodinger) 
semigroup - driven dynamics. It turns out that the semigroup modeling provides much stronger 
confining properties than the standard Langevin one. Since contractive semigroups set a link be- 
tween Levy flights and fractional (pseudo-differential) Hamiltonian systems, we can use the latter 
to control the long - time asymptotics of the pertinent pdfs. To do so, we need to impose suitable 
restrictions upon the Hamiltonian and its potential. That provides verifiable criteria for an invariant 
pdf to be actually an asymptotic pdf of the semigroup-driven jump-type process. For computational 
and visualization purposes our observations are exemplified for the Cauchy driver and its response 
to external polynomial potentials (referring to Levy oscillators), with respect to both dynamical 
mechanisms. 
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I. INTRODUCTION 

In recent years there has been growing interest in ran- 
dom walks, extending from various fields of physics to 
chemistry, biology and financial mathematics. The clas- 
sical concept of Brownian motion has become paradig- 
matic in the whole theory of stochastic processes, see 
e.g. Ref. The probability density function (pdf) of 
a homogeneous Brownian motion solves a Fokker-Planck 
equation and has an important intrinsic property: the 
diffusing particle pdf, that is initially concentrated at a 
point, with the flow of time takes the Gaussian form, 
whose width grows in time as t 1 / 2 . This kind of dif- 
fusion processes was called the normal diffusion. We 
leave aside a broad field of anomalous diffusions, where 
< X 2 (t) >~ t a with < a < 2, and focus on non- 
Gaussian, jump-type stochastic processes, whose pdfs be- 
long to a class of Levy - stable distributions. They have 
long tails and their second (and higher) moments are 
nonexistent. The existence of the first moment is granted 
only for a suitable subclass Q. 

In the standard theory of diffusion- type processes, a 
celebrated method to solve transport problems with re- 
spect to the pdf and specifically transition probability 
density, is to transform the involved Fokker-Planck equa- 
tion into its parabolic (Hermitian) counterpart, Right 
at this point a dynamical (Schrodinger) semigroup notion 
naturally appears and its obvious link with a Hamilto- 
nian dynamical system. Indeed, one encounters here the 
Schrodinger-type equation (there is no imaginary unit 
before time derivative and physical dimensions of the 
Hamiltonian are re-scaled accordingly), often called the 
generalized diffusion equation, dt'F = H^f. Here, H is 
interpreted as a Hamiltonian operator and — H stands 
for a generator of the dynamical semigroup exp(— tH). 

In case of Brownian motion and normal diffusion pro- 
cess, the Fokker - Planck (Langevin - driven) and semi- 
group dynamics refer to the same random process. A use- 
fulness of the the semigroup picture lies in the possibility 
of eigenfunction expansions of H which allows to deduce 



explicit formulas for transition (semigroup) kernels and 
transition pdfs, [3]. As a byproduct of the above proce- 
dure one solves the eigenvalue problem for the Hamilto- 
nian operator and identifies its ground state as a square 
root of an invariant pdf of the stochastic diffusion pro- 
cess. The latter pdf is approached by the process at large 
times. 

The situation is somewhat different in the case of 
non-Gaussian jump-type processes. Namely, the (Levy 
- Schrodinger) semigroup dynamics differs from the 
Langevin - driven fractional Fokker-Planck evolution, see 
e.g. 0, H, @|. Nonetheless a common asymptotic invari- 
ant pdf may be attributed to both dynamical scenarios 
4]. Since Langevin-driven dynamics does not admit in- 
variant pdfs in the Gibbs form (e.g. the force potential 
does not appear in the exponential form of p*), the whole 
class of pdfs, originally employed in the study of topo- 
logically induced super-diffusions [5|, [fjl , might seem to be 
excluded by the formalism of Ref. Q. This is not the 
case. 

In the present paper we focus on extending the range 
of validity of the reverse engineering (targeted stochas- 
ticity) problem of Ref. [?} to so-called topological Levy 
processes (topologically induced super-diffusions), occur- 
ring in systems with topological complexity like folded 
polymers and complex networks. To be more specific, 
the original reverse engineering problem formulates as 
follows: given an invariant pdf p*(x), design a stochastic 
Langevin-driven jump-type process for which the prese- 
lected density may be an asymptotic target. The basic 
reconstruction goal is to deduce the drift function of the 
process. 

In the previous paper, Q, we have recast the re- 
verse engineering problem so that the original task 
was supplemented by one more reconstruction step. 
Namely, we have addressed the existence issue of 
the Levy - Schrodinger semigroup potential V(x) = 

—A (| A\ tJ, ^ 2 p 1 J 2 )(x) /p 1 J 2 (x) (see below for detailed expla- 
nation), given the very same (as for the Langevin pro- 
cess) invariant pdf p*(x), that is non-Gibbsian by con- 
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struction. 

Presently, we relax the previous (common pdf) con- 
straint and address a fully fledged reconstruction prob- 
lem for the semigroup dynamics: given an invariant pdf, 
identify the semigroup-driven Levy process for which the 
prescribed invariant pdf /?* may stand for an asymp- 
totic one. Then, Gibbsian densities appear to be admis- 
sible and a class of jump-type processes, that respond 
to environmental inhomogeneities, becomes largely ex- 
tended. The corresponding jump- type processes are iden- 
tified by us as topological processes due to their links with 
topologically- induced super-diffusions, d, [1]. 

At this point lest us stress that it has never been settled 
that the invariant pdfs of a topological process actually 
are the proper asymptotic ones, e.g. can be reached in 
the large time asymptotics irrespective from a particular 
choice of initial data. To this end one must resort to the 
contractive semigroup notion. One of the aims of our 
present discussion is to carefully check this point. 

The expected asymptotic behavior may not persist 
for an unrestricted initial pdf choice. The signature, of 
whether such behavior is allowed or prohibited by the 
semigroup dynamics, is encoded in the functional form 
of a semigroup potential V(x). The latter needs to be 
reconstructed from the target pdf p*(x), by means of 
the above generalized reverse engineering problem and 
respect a number of restrictions. Minimal requirements 
upon the associated pseudo-differential Hamiltonian and 
its potential V(x) were set in Ref. [To[, where an ex- 
plicit construction of Cauchy semigroups has been car- 
ried out. First explicit examples of appropriate poten- 
tials were found in Q. It is the contractive semigroup 
dynamics that guarantees a proper asymptotic behavior 
of inferred time-dependent pdfs, c.f. also [8(] for a more 
advanced exposition of that issue. 

II. LEVY SEMIGROUPS IN A RANDOM 
MOTION 

A. Brownian pre-requisites 

If we have a one-dimensional Smoluchowski diffusion 
process [1[ with an initial pdf po(x), then its time 
evolution is determined by the Fokker-Planck equation 
dtp = DAp — V (b ■ p) where D is a diffusion coeffi- 
cient and the time - independent drift b(x) — f(x)/m(3 = 
— (l/m/3)VV(x) is induced by an external (conservative, 
Newtonian) force field f(x) — — VV(x). We adopt a 
standard form D = k B T jmfi of the diffusion coefficient, 
where m and (3 are, respectively, a mass and a reciprocal 
relaxation time of a particle. 

Following a standard procedure [|| we may identically 
rewrite the Fokker-Planck equation in terms of an asso- 
ciated Hcrmitian (Schrodinger-type) problem by moans 
of a redefinition 

p(x,t) = V(x,t)pl /2 (x) (1) 



that takes the Fokker-Plack equation into a parabolic 
one, often called a generalized diffusion equation: 

df* = DA# - V* (2) 

for a positive function ty(x,t). The auxiliary poten- 
tial V derives from a compatibility condition V(x) = 
DAp^J 2 / ' pV 2 \ whose equivalent form reads V(x) = 
{\/2)[&/{2D) + Vb]. 

If the {\/2mD rescaled) Schrodinger-type Hamiltonian 
H = —DA + V is a self-adjoint operator in a suitable 
Hilbert space, then one arrives at a dynamical semigroup 
exp(— tH), We note here, that the Schrodinger semi- 
group (parabolic) reformulation of the Fokker - Planck 
equation is merely another mathematical "face" of the 
diffusion process, the operator H is just one more form 
of the Fokker-Planck operator [f| . The semigroup is con- 
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tractive, hence asymptotically ^(x, t)\t— ,00 ~ > P* {%). 
Accordingly, p{x,t)\ t ^ SXl — > p*{x). 

We note that for V = V(x) bounded from below, 
the integral kernel k(y,s,x,i) — {exp[— (t — s)H]}(y,x), 
s < t.of the dynamical semigroup exp(— tH), is pos- 
itive Q. The semigroup dynamics reads: \£(x, t) — 
j 5 , (y, s) k(y, s, x, t) dy so that for all < s < t 

1/2 / 
p(x,t) = pj (x)V(x,t)= p(y,s,x,t)p(y,s)dy, (3) 

where 

p 1 ^ 2 (x) 

p(y, s, x, t) = k(y, s, x, t)^-^ — (4) 

p* (y) 

is the transition probability density of the pertinent 
Markov process. Its unique asymptotic invariant pdf is 
P*{x). 

For the familiar Ornstein-Uhlenbeck version of the 
Smoluchowski process, the drift is a linear function of 
x, e.g. b(x) = —72;, 7 = K/m/3, k > 0. The Fokker - 
Planck equation dtp — DAp + ^(x p) supports an in- 
variant density 




where V(x) = k^- and F* = -k B T\rL{2irk B T / n) 1 ' 2 . 

The associated generalized heat equation involves H = 

2 2 

—DA + V with V(x) = 2jg ^ which is is a typical 

1/2 

confining potential. Accordingly, HpJ = 0. The lowest 

eigenvalue of this positively defined operator identifies 

1/2 

pj as its ground state function. Eq. (fj| rewrites as 
p(y, s, x, t) = k{y, s, x, t) exp[V(y) - V{x)]/2k B T . (5) 

B. Free Levy-Schrodinger Hamiltonian 

To consider the properties of a free (without external 
potentials) Levy-Schrodinger semigroup, we employ the 
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rescaled Hamiltonians rather than semigroup generators 
that have an opposite sign. The pertinent Hamiltoni- 
ans have the form H = F(j>), where p — — iV stands 
for the momentum operator (up to the scaled away H 
or 2mD factor), and for — oo < k < +00, the func- 
tion F = F(k) is real valued and bounded from be- 
low. The action of exp(-tH) can be given by means 
of an integral kernel k t = k(x — y,t) — fc(y, 0, x, t) where 
fct ( z ' *) = ^fc f-00 exp[-tF(p) + ipz]dp. 

Our further discussion is limited to non - Gaussian 
random variables whose pdfs are centered and symmetric, 
e.g. to a subclass of stable distributions characterized by 



F(p) = Ab|' 



A|A|^ 2 . 



(G) 



Here < p < 2 and A > stands for an intensity parame- 
ter of Levy process. To account for the interval < p < 2 
boundaries, one should rather employ (—A)''/ 2 instead of 
|A|^/ 2 ; —A is a positive operator. 

The pseudo-differential Hamiltonian by construc- 
tion is positive and self-adjoint on a properly tailored do- 
main. A sufficient and necessary condition for both these 
properties to hold true is that the pdf of the Levy process 
is symmetric, see Ref. [l3j . The corresponding contrac- 
tive semigroup admits an analytic continuation in time 
leading to Levy-Schrodinger equations and fractional 
quantum mechanics [IJ, [IH, [Hj]. The associated jump- 
type dynamics is interpreted in terms of Levy flights. The 
pseudo-differential Fokker-Planck equation, which corre- 
sponds to the free fractional Hamiltonian and the 
fractional semigroup exp(— tH^) = exp(— A| A| M / 2 ), reads 



d tP = ~\\Ar/ 2 P , 



(7) 



to be compared with the ordinary heat equation dtp = 
+DAp. In particular F(p) = X\p\ refers to Cauchy pro- 
cess. 

The action of the pseudo-differential operator | A\^ 2 
on a function can be expressed by the formula [lol [l2| 

{\A\^f){x) = - I [f(x + y)- f(x)]^(dy), (8) 



where v^(dy) is a corresponding Levy measure (see, e.g. 
0) and the integral in Eq. is understood in a sense 
of its Cauchy principal value. Changing the integration 
variable y to z = x+y and employing a definition of Riesz 
fractional derivative of the p-th order, [l^|, we arrive at 



r( M +l)sin(7r M /2) ff(z)-f(x) 



m^mx) = — j f^f^z 

(9) 

with (|A|^/ 2 /)(a;) = -d»f{x)/d\xf. The case of p = 1 
refers to the Cauchy driver (e.g. noise). We note a 
systematic sign difference between our notation for a 
pseudo-differential operator | A| p / 2 and that based on the 
fractional derivative notion, like e.g. A M / 2 = /d\x\^ of 
Refs. HQ. 



C. Response to external potentials 

Consider now the Levy-Schrodinger Hamiltonian with 
external potential 



(10) 



Suitable properties of V need to be assumed, so that 
— Hp is a legitimate generator of a dynamical semigroup 
exp(— tH^), see e.g. Ref. [Toj . 

Looking for stationary solutions of the equation dflf = 

H^, we realize that if a square root of a positive in- 

1/2 

variant pdf ^ ~ pj is asymptotically to come out, then 

fractional Sturm-Liouville operator should be used to de- 
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rive an explicit form of pj for a given V. In the opposite 
situation, when p*{x) is a priori prescribed, we can de- 
termine V through a compatibility condition: 



V = -A 



|A|"/ 2 pi 



1/2 



P 



1/2 



(11) 



The main point here is that we do not have too much 
freedom in pre-selecting a functional form of p*, as suit- 
able conditions need to be respected by the inferred aux- 
iliary potential V, to yield a contractive semigroup dy- 
namics. This leads to a conclusion [Io| that only under 
the contractive semigroup premises, an invariant pdf of 
a jump-type process may actually become its asymptotic 
target. The detailed discussion of this issue for Cauchy 
semigroups can be found in Ref. 10 . 

To derive the pseudo-differential equation governing 
the behavior of a system in an external potential, we 
rewrite the pdf of the semigroup-driven stochastic pro- 
cess in the form ([T]) (see also [4]). Any strictly positive 
function (here we consider only functions that vanish for 
large x) can be rewritten in an exponential form. Hence, 
by adopting the notation p*(x) = exp[2<I>(a;)] and ac- 
counting for (jlip we arrive at a continuity equation with 
an explicit fractional input 



d tP = -A(exp$)|A| /i/2 [exp(-$)p] + Vp. 



(12) 



The definition (jTTJ) suggests that any pdf of the form 
p*(x) = exp[2$(x)] is a proper candidate for a stationary 
solution of the Eq. ([12]) . However things are not that 
simple. It is only the semigroup dynamics generated by 
(Tll"j) and (fl"2"|) that may guarantee a consistent temporal 
approach towards an asymptotic invariant density of the 
stochastic process in question. Right at this point, an 
issue of restrictions upon an effective potential V of Eq. 
(dl]), dEl, enters the stage. 

It is instructive to mention that for Levy flights in ex- 
ternal force fields, the (somewhat left aside) Langevin 
approach is known to yield [Tt} a continuity (e.g. frac- 
tional Fokker-Planck) equation in a very different form 



dtp = -V 



H S p\-\\M" 2 , 

mp 



(13) 
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Even if we know [i[ , that an asymptotic invariant density 
ofEq. (|12p may coincide with that for Eg. (fT5)) . these two 
transport equations refer to different temporal patterns 
of behavior. 

We note that, in contrast to the semigroup modeling, 
the Langevin scenario for Levy flights in confining po- 
tentials has received ample attention in the literature, 
see @, [U El [H H3 to cite a few. 

D. Topologically induced super-diffusions 

It is of some interest to invoke an independent (so- 
called topological, see above) approach of Refs. [1, @ 
where one modifies jumping rates by suitable local factors 
to arrive at a response mechanism that is characteristic 
of the previously outlined semigroup dynamics. Namely, 
in view of (9), the free transport equation (7) can be 
re-written as a master equation: 

d t p{x) = J[w(x\z)p(z) - w(z\x)p(x)]v fl (dz) . (14) 

The jump rate is an even function, w(x\z) = w(x\z). 
However, if we replace the jump rate 



w(x\y) ~ l/\x~y\ 1+ » 
(c.f. Eq. ©) by the expression 

exp[$(x) - $(y)] 



w<t>( x \y) 



\x - y\ 



(15) 



(16) 



J 



and account for the fact that w c f ) (x\z) ^ w c f,(z\x), the 
corresponding transport equation takes the form: 



(1/X)d t p = |A|£ /2 / = - exp($) |Ar/ 2 [exp(-*)/ 
+pexp(-$)|A|' i/2 exp($). 



(17) 



Whatever potential <&(x) has been chosen (up to a nor- 
malization factor), then formally p*(x) = exp(2$(x)) is 
a stationary solution of Eq. (|I7[) . Moreover, one readily 
verifies that Eq. (|I7p is identical with the semigroup- 
induced Eq. (fT2]) . 

Accordingly, if for a pre-determined = exp(2$), 
there exists the semigroup potential V, Eq. (fTTj) . then 
Eq. (fTT)) defines the dynamics that belongs to the pre- 
viously outlined semigroup framework. This entails a 
direct verification of whether the stationary density 
may really be interpreted as an asymptotic target of the 
pertinent fractional transport equation. 

Stationary pdfs for the topologically-induced dynam- 
ics were demanded to occur in a Gibbsian form exp(2$), 
see 0, Q for a possible phenomenological background 
for this assumption. Therefore, we redefine $ as fol- 
lows. Rewriting the stationary pdf as p*{x) = 
(1/Z) exp(— V*(x)/kBT) (normalization factor Z stands 
for a partition function), we recover a function V*(x) — 
— ksT \n(Z p*(x)) that receives an interpretation of the 
external potential. With these re- definitions, the previ- 
ous equation (|17p takes the customary form employed in 
the discussion of topologically induced super-diffusions: 



dtp = -exp(-/tK/2) |A|^ 2 exp(KK/2)p + pexp( K V;/2)|A|' J / 2 exp(-KK/2), k = l/k B T. 

I 



(18) 



III. RESPONSE OF CAUCHY NOISE TO 
POLYNOMIAL POTENTIALS 

A. Ornstein - Uhlenbeck - Cauchy process 

In case of the Ornstein - Uhlenbeck - Cauchy (OUC) 
process, the drift is given by b(x) — —jx, and an asymp- 
totic invariant pdf associated with 



d t p = -\\V\p + V[frx)p] 



reads: 



p*(x) 



a 



7r(cr 2 + a; 2 ) ' 



(19) 



(20) 



where a = A/7, c.f. Eq. (9) in Ref. [1JJ. 

Note that a characteristic function of this density reads 
F(p) = ~ a \p\ an d gives account for a non - thermal fluc- 
tuation - dissipation balance. The modified noise inten- 
sity parameter a is a ratio of an intensity parameter A of 
the Cauchy noise and of the friction coefficient 7. 



The invariant density of OUC process (|20[) generates 
(with the help of Eq. (fTTj) ') the following Cauchy semi- 
group potential V (we set p. = 1) 



V(x) = - 

IT 



'a 



x 



x 1 
-In 

a s/a — x 



(21) 



The potential (|2f|) had been analyzed in Ref. [4]. It is 
clear that V(x) is bounded both from below and above, 
well fitting to the general mathematical construction of 
(semigroup-driven or topological) Cauchy processes in 
external potentials, [IgT |. 

Since the OUC pdf ([20]) has no variance, in Fig. [JJ 
we visualize the temporal evolution of OUC process with 
initial data localized at x = (p(x, t = 0) = 5(x)) in two 
motion scenarios (i.e. Langevin and semigroup driven) 
by comparing the width of the OUC "bell" at its half- 
maximum at a number of consecutive instants of time. It 
is seen that Langevin dynamics sets at equilibrium faster 
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FIG. 1: Temporal behavior of the half-maximum width (HW): 
for the OUC process in Langevin-driven and semigroup-driven 
(topological) processes. Motions begin from common initial 
data p(x,t = 0) = 8(x) and end up at a common pdf (|2Up for 

(7=1. 



FIG. 2: Time evolution of Langevin-driven pdf pL{x,t) be- 
ginning from the initial data pL(x,t = 0) = 5(x + 1) and 
ending at the pdf (|20|) (shown as "asymptote" in the figure) 
for a — 1. Figures near curves correspond to t values. 



B. Cauchy driver: Polynomial drifts and the 
semigroup evolution 



than the semigroup-induced dynamics. 

Further check of the OUC dynamics is to obtain (nu- 
merically [2l|) the time evolution p(x,t) of initial data, 
localized in another point. Since quadratic potential (cor- 
responding to OUC process) has a minimum at x = 0, 
we investigate the temporal evolution for the initial data 
that are shifted from the point x — 0. Namely, we con- 
sider two cases - "left" pl(x, t = 0) = 5(x+l) and "right" 
p R (x,t = 0) = S(x - 1) ones. As the "left" and "right" 
cases are symmetric with respect to y axis, in Fig. [5] we 
report the "left" case only. It is seen that final stage of 
the evolution is still the invariant pdf ([20]) . This means 
that fictitious particle representing our process, "rolls 
down" to potential minimum at x = 0. We note that this 
" rolling down" occurs slower then the evolution with ini- 
tial data, localized at x = 0. The topologically- induced 
dynamics of pL t R{x,t = 0) is qualitatively the same, 
except for the just mentioned slow-down, as Langevin- 
induced one (seeRef. (|) and we do not show it here. 

A pre-selection of the OUC p» , Eq. ([2"U|) , proved to pro- 
vide (through a numerical process reconstruction) a con- 
sistent asymptotic target pdf for the Cauchy semigroup- 
driven dynamics. At this point we may justifiably ask 
about limitations upon freedom of choice that is admit- 
ted in such invariant pdf pre-selection procedure. 

Remark 1: As a byproduct of the discussion we 

have established a pseudo-differential Hamiltonian sys- 

l li 

tern, whose ground state equals the square root pj (x) 
of the Cauchy pdf (H|). 



The OUC case corresponds to a linear drift function. A 
number of polynomial drift functions has been discussed 
in the literature with a focus on confining features of 
various external forces on Levy flights. In each case a 
corresponding asymptotic invariant pdf has been found, 
albeit with a restriction (in view of limited computational 
facilities) to the Cauchy driver. For clarity of subsequent 
discussion below we provide some examples. 

The quadratic Cauchy density p*(x) = 2/n(l + x 2 ) 2 
stands for an asymptotic target of the Langevin - driven 
process with the drift b(x) = (~jx/8)(x 2 + 3) 4]. For the 
above quadratic Cauchy pdf, the associated semigroup 
evolution is defined, by means of the potential function 
V(x) = \{x 2 - l)/(x 2 + 1), obtained from This 
potential is also bounded from below and above and con- 
sequently, [l(| , defines a contractive semigroup generator 
and hence semigroup dynamics. 

The bimodal pdf pj{x) = (3 3 /tt(x 4 - (3 2 x 2 + /3 4 ) can 
be derived from the drift function b{x) = —jx 3 /f3 3 . The 
associated semigroup dynamics is correctly defined by 
means of a potential function V, obtained from Eq. 
numerically. Its shape is reported in Fig. [3J It is seen, 
that the potential is bounded from below and above, 
again perfectly fitting to the general theory of [l(| • 



C. Cauchy driver: Gibbsian versus non-Gibbsian 
asymptotics 

The Langevin-driven dynamics with a given drift b ~ 
— VV, where — VU stands for a conservative force acting 
upon a particle in the course of its random motion, for 
non-Gaussian driver (like e.g. Cauchy one) does not ad- 
mit an asymptotic invariant density in the Gibbsian form 
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FIG. 3: The coordinate dependence of V(x) associated with 
the bimodal pdf, for different values of f3. 



p*~exp{-V/k B T), [7|. 

As we have established above, see also the 
Langevin-driven and semigroup-driven jump-type pro- 
cesses, with the Cauchy driver, may share common 
asymptotic stationary pdfs, that are obviously non- 
Gibbsian. 

On the other hand, as mentioned in Section II. D, 
one may suspect that asymptotic invariant pdfs of a 
semigroup-driven process (in other words, topologically- 
induced one) may not necessarily have fat (e.g. non- 
Gaussian) tails, due to extremely strong confinement 
of admissible jumps, imposed by the "potential land- 
scape" of an inhomogeneous medium. That amounts 
to assuming that a Gibbs-type asymptotic pdf p* = 
(l/Z) exp(— V^/ksT) may be employed in the construc- 
tion of the topologically- induced dynamics (fPSj) . 

A thermalization mechanism, under which the equilib- 
rium conditions could have been achieved for such Gibb- 
sian (the non-Gaussian mechanism is excluded 0] and 
has not been considered in Refs. @, 0) is actually un- 
clear. The major focus in the literature was upon a suit- 
able "potential landscape" (potential profile), such that 
local modifications (|16p of the jump rates would drive the 
random motion towards a Gibbsian equilibrium. 

We shall follow the " potential landscape" intuition and 
refer to explicit dynamics simulations reported in Ref. [fj 
for the double well potential V*{x) = — x 4 — 2x 2 + 
1. Dimensional units are scaled away and, to facilitate 
comparison, the notation $ refers presently to that of 
Ref. @), see Fig. 3 therein. As a supplementary test, 
we consider also $ = V*(x) = x 2 , leading to Gibbsian 
asymptotic pdf in the Gaussian form. 

Our next step was to evaluate (numerically) the semi- 
group potential V, after literally substituting a concrete 
Gibbsian asymptotic pdf p* to Eq. (fTTj) . The outcome is 
reported in Fig. 2J 

The semigroup potentials V depicted in Fig. 0] are 



FIG. 4: The coordinate dependence of the semigroup po- 
tential V(x) (curves 1 and 2), corresponding to V,(x) = 
x 4 — 2x' 2 + 1 (curve 3) and V* (x) = x 2 (curve 4), respectively. 
Curves 3 and 4 are shown for a comparison with, strikingly 
similar in shape, semigroup potential curves 1 and 2 



admitted by the general theory of Ref. [lCj (see below for 
details), so granting the operator exp(—tH) (H — — | V| — 
V), the contractive semigroup status. This proves that 
an invariant (stationary) density of Eq. (|19[) . with the 
bimodal or harmonic exponent, actually is an asymptotic 
invariant pdf of the topological process. 

The topological confining mechanism appears to be 
much stronger than confining mechanism based on the 
Langevin modelling. Namely, if both mechanisms share 
the same potential $>(x) — V*(x), then the asymptotic 
pdf in Langevin scenario has no more then a finite num- 
ber of moments. At the same time, the corresponding 
asymptotic pdf of the topological process, being in the 
form cx exp(— $), under the Gibbsian premises may in 
principle admit all moments. To the contrary, if we im- 
pose for both mechanisms to have the same asymptotic 
pdf, the above statement is invalid. 

For completeness, we now recollect the above men- 
tioned requirements upon V that need to be observed in 
the presence of Cauchy driver, [l(|- Namely, the poten- 
tial should allow to be made positive (by merely adding 
a constant), should be locally bounded and needs to be 
measurable (i.e. can be approximated with arbitrary pre- 
cision by step functions sequences). The Cauchy gener- 
ator plus a potential with such properties is known to 
determine uniquely [l(| an associated Markov process of 
the jump-type and its step functions approximation. The 
limiting behavior of the pertinent step process, as the 
step size is going to zero, remains under control. 

Remark 2: Would we have followed the standard 
Langevin modeling for the Cauchy driver, with the exter- 
nal force potential V* (x) = x 4 — 2x 2 + 1 and the resultant 
drift — VT4 = b, an invariant pdf of the corresponding 
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FIG. 5: The coordinate dependence of the semigroup poten- 
tial V(x) derived from the non-Gibbsian pdf of Remark 2, 
Section III.C, along with the pdf J22J. 

fractional Fokker-Planck equation would have the form: 
_ 2a(a 2 + b 2 ) 1 

P*W ~ ~ ( a 2 +6 2)2 +2 ( a 2_ 6 2) x 2 +:E 4 

(22) 

with a ~ 0.118366 and b ~ 1.0208. Here, a and b are, 
respectively, real and imaginary parts of complex roots 
of the cubic equation z 3 + z — 1/4 = 0. It is easy to 
show both analytically and numerically that the above 
p* is properly normalized. By using the formula (jllj) we 
may associate with the non-Gibbsian pdf a semigroup 
potential. Its properties prove that we deal with a con- 
tractive semigroup dynamics and a common asymptotic 
pdf for both Langevin and semigroup motion scenarios. 
The outcome of this discussion is depicted in Fig 5. 
Remark 3: Standard arguments , c.f. Section II. A, con- 
vince us that the Gibbsian form p* ~ exp(— V^/ksT) 
of an asymptotic density would have been recovered in 
the presence of the Wiener (Brownian) noise, given the 
drift function b ~ — VVi with V„{x) = x 4 — 2x 2 + 1 or 

V*(S) = X 2 . 

D. Cauchy oscillator 

Our preceding discussion has actually been related to 
the generalization of the reverse engineering problem of 
Ref. [7J: (i) with an invariant pdf p*(x) in hands, de- 
rive b(x) for the associated Langevin equation to recon- 
struct the Langevin-drien dynamics of the pdf, (ii) from 
the same invariant pdf p*(x), deduce the topologically - 
driven potential V(x) for the related Cauchy semigroup 
dynamics (thus granting that pdf an asymptotic pdf sta- 
tus). 

In the present section we invert the reasoning and ef- 
fectively follow the direct engineering route: having V(x), 



we employ the semigroup dynamics principles to infer an 
asymptotic target p*(x). This construction will be sup- 
plemented by identifying, if any, the drift function b(x) 
for the associated Langevin process that gives rise to the 
same asymptotic target. 

To this end it seems natural to resort to simplest pos- 
sible potential functions. Our choice (motivated by a 
simplicity of involved Fourier transforms and an imme- 
diate comparison with the Gaussian OU process of Sec- 
tion II. A) is the familiar harmonic oscillator potential 
V(x) = § x 2 — V , constant Vo is left unspecified. 

Our major object of interest is thus a pseudo- 
differential Hamiltonian (this is a very special, massless 
version, of the well known Hamiltonian for the relativistic 
harmonic oscillator problem, [l2l [l3l]): 

H 1/2 = A|V| + (|x 2 - Vo) . (23) 

That in turn is a Cauchy analog of the familiar harmonic 
oscillator Hamiltonian of Section II. A: 

*— "^Ossr-?)- < 24 » 

Since the quadratic function is admissible 10] as a semi- 
group potential, we need not to bother about an asymp- 
totic approach towards an invariant density, but follow a 
direct reconstruction route: given V of Eq. (f2"4"]l . deduce 
an invariant pdf . To this end we turn back to Eq. (jll[) 
with p — 1 and consider a pseudo-differential equation 

(f^V ) P y 2 = -A|V|pi/ 2 (25) 
to solve it with respect to 

~ 1/2 

We denote f(p) the Fourier transform of f — pj (x). 
Accordingly, Eq. takes the following form: 

-~A p f + j\p\f = V f (26) 

which, up to constants adjustments, can be recognized as 
the eigenvalue problem for the Schrodinger operator with 
the linear confinin g (m odulus of the argument) potential, 
[22| . see also [23l |24J. l25j . Albeit with respect to mo- 
mentum - space (here, wave-vector) variables, i.e. with 
A p = d 2 /dp 2 replacing the conventional spatial Lapla- 
cian. 

By changing an independent variable p to k — (p — 
next denoting ijj(k) — f(p) with the identifications 
(j = V0/7 and C = (K/27) 1 / 3 , we may rewrite the above 
eigenvalue problem (with Vo standing for an eigenvalue) 
in the form of the following ordinary differential equation 

^ = mm, (27) 

whose solutions can be represented in terms of Airy func- 
tions. A brief resume of how to deduce the eigenfunctions 
and eigenvalues of the original problem (|26p , is relegated 
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x 



FIG. 6: Fourier image t[>o(k) of the function ij)o(x) along with 
this function 



to Appendix. A unique normalized ground state function 
of the problem (f2"rJ)) , (J27J) is composed of two Airy pieces 
that are glued together at the first zero yo of the Airy 
function derivative: 



ip (k)=A 



An 



Ai(-j/ + k), k>0 
Ai(-y ~k), k< 0, 



y w 1.01879297. (28) 



The function Q28p is a square root of the pdf in momen- 
tum space. We reproduce its shape in Fig. [6j 

To transform the ground state solution back to coordi- 
nate space, we evaluate the inverse Fourier transforma- 
tion of the ground state solution (|2"8")l . see Appendix C 
for details. This yields the following real ground state 
wave function f(x) — > iI>q{x) 



A 



ip (x) = — I Ai(t) cos x(t + y )dt = pl /2 (x), (29) 
n J-y 

which determines an invariant pdf p*(x) of the direct 
engineering problem (|25p as follows: 




FIG. 7: Normalized invariant pdf J30} (full line) for the 
quadratic semigroup potential. The Gaussian function, cen- 
tered at x = and with the same variance a 2 = 0.339598 
is shown for comparison. Panel (a) shows functions in lin- 
ear scale, while panel (b) shows them in logarithmic scale to 
better visualize their different behavior. 



p*{x) 



Ac 



Ai(i) cosx(t + yo)dt 



(30) 



Ao 

7T 



dt / dtiM(t)M(ti) x 
-yo J —yo 
x cosx(t + y ) cosx(h + y ). 



a 2 = 0.339598, f(x) 



' . It is shown by a 



dashed line. Fig. [7Jd shows both above functions in the 
log scale. It is clearly seen, that Airy- induced p*(x) de- 
cays at infinities slower then Gaussian. 



Even with an exact analytic formula for the normalized 
p*(x)dx = 1 function p*(x) (|30|) in hands, a better 
insight into its properties is achieved only by means of 
numerical methods. We depict both x and p-space ver- 
sions of "00 m Fig. 6. 

In Fig. 7, the resultant invariant density p*(x) is rep- 
resented by a full line. For comparison, we have de- 
picted the Gaussian function with the same variance 



Let us finally note, that the potential $(x) may be 
recovered from the ground state function (|30|) as the 
asymptotic pdf for topologically driven (semigroup) pro- 
cess has the form exp(— $). Namely, $(x) oc — lnp*(x), 
where p*(x) is given by Eq. (|3TJ)) . The inversion of the 
y axis of Fig. [TJd shows this potential (compared with 
harmonic one y = x 2 ). 
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FIG. 8: Langevin - type drift b(x) (curve 1) and its (force) 
potential V(x) (curve 2), that give rise to an invariant density 
(30). 



E. Reverse engineering for the Cauchy oscillator 
ground state pdf 

For a given p* the definition of a drift function b(x) 
(we put either A = 1 or define b — > b/X) is: 



b(x) 



p*(x) 



\V\p*(x))dx : 



(31) 



(x) 



dx 



p*(x + y) - p*(x) 



dy. 



Inserting p*(x), Eq. ([5T3|) . we get 



b(x) = - 



C yo M(t) S mx(t- 



y )dt 



I- yo Ai W cos x(t + y )dt' 



(32) 



The final formula for b(x), (f3"2"|) , together with that 
for the corresponding Langevin force potential V(x) = 
— J b(x)dx look a bit clumsy. Therefore it is appropri- 
ate to reiterate to numerics, see Fig. El The plot of 
b(x) is reported in FigJS] along with potential function 
V{x) = -fb(x)dx. 

We were unable to determine the asymptotic of b(x) at 
spatial infinities. It is possible, however, to expand Eq. 
(f3"2"| in power series at small x. It turns out, that these 
truncated series describe the function b(x) (and corre- 
spondingly V(x)) surprisingly well. 

To obtain this approximation, we expand both nu- 
merator and denominator of Eq. (|32p in power se- 
ries to obtain Ai(t) smx(t + y )dt = 0.824278x - 

0.503237x 3 + 0.205648a; 5 -0.0650381x 7 +... for numerator 
and f™ g Ai(i) cos x(t + y )dt = 0.809073 - 0.687715a; 2 + 

0.334243a; 4 - . 1 1 8792a; 6 + 0.0339885 x 8 - . . . for denomina- 
tor. The series can be easily continued for larger amount 



of terms. Now, the ratio of these series can be also ex- 
pressed in the form of the series. Here, we reproduce only 
the truncated (up to x 5 ) version of the series 

b(x) w -1.01879a:- 0.243986a: 3 - 0.04056x 5 , 
V(x) w 0.509395a; 2 + 0.0609965a; 4 + 0.00676x 6 .(33) 



The approximations (|33[) can be used to calculate nu- 
merically the Langevin-type dynamics with the invariant 
density f30]). 



IV. CONCLUSIONS 

The Levy - Schrodingcr semigroup modeling sets a 
(mathematically rigourous) link between Levy flights and 
pscudo - differential Hamiltonian systems. It is well 
known that the pdf of a free Levy flight has so-called 
heavy tails and consequently does not have second and 
higher moments. Properly tailored external inhomo- 
geneities are capable of "taming" a Levy flight so that 
the resultant pdf admits higher moments which, in turn, 
makes possible to use such functions for the description 
of real physical (and other, like biological or economic) 
systems. 

We have encoded the overall impact of inhomogeneities 
in the semigroup potential notion, that may be inter- 
preted as an external potential in the affiliated pseudo- 
differential Hamiltonian operator. However, not any con- 
ceivable semigroup potential and thus, not any conceiv- 
able inhomogeneity, can make the corresponding jump- 
type process a mathematically well-behaved construc- 
tion, where the Markovian dynamics entails an approach 
towards a unique stationary pdf. 

In the present paper, we have investigated the behav- 
ior of Levy oscillators in different external confining po- 
tentials. Our analysis shows that to control the long - 
time asymptotics of the above Levy oscillator pdfs, suit- 
able restrictions upon the Hamiltonian and its (semi- 
group or topologically-induced) potential V{x) need to 
be observed, [l(J. 

Namely, V(x) should allow to be made positive (this is 
achieved by simple vertical shift of the entire function), 
should be locally bounded and needs to be measurable, 
i.e. should have a possibility to be approximated with ar- 
bitrary precision by step functions sequences. The fulfil- 
ment of these requirements provides verifiable criteria for 
an invariant pdf to be actually a time - asymptotic pdf of 
a semigroup (equivalently, topologically)-driven process. 

A technical advantage of a semigroup formalism is the 
possibility of eigenfunction expansions for H which allows 
to deduce explicit formulas for transition pdfs, As a 
byproduct of the above procedure we have completely 
solved the eigenvalue problem for the Hamiltonian oper- 
ator of so-called Cauchy oscillator in terms of Airy func- 
tions. The ground state wave function of such oscillator 
has been obtained analytically. Its square defines the in- 
variant pdf of Cauchy oscillator process. The latter pdf 
is approached by the jump-type process at large times. 
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We have extended the targeted stochasticity problem 
of Ref. 0] to the above semigroup-driven (topological) 
Levy processes, which are widely used in literature to 
model various systems, like polymers, glasses and com- 
plex networks. Our departure point was as follows: hav- 
ing an invariant pdf p* (x) , recover not only the Langevin 
drift b(x) and potential V(x) = — J b(x)dx, but also the 
potential V(x) of the corresponding topological (semi- 
group) Levy process, being attributed to the same in- 
variant pdf. 

Furthermore, we have relaxed a common pdf require- 
ment and have reformulated the targeted stochasticity 
problem as a task of reproducing a suitable contractive 
semigroup, given an invariant pdf, with the Levy (specif- 
ically, Cauchy) driver in action. We have shown, that 
the semigroup modeling provides much stronger confin- 
ing properties than the standard Langevin one, such that 
the resultant asymptotic pdf may have all moments. 

To be more specific, if both above approaches involve 
(albeit differently) the same conservative force potential 
$ = V*(x), then the asymptotic pdf in the Langevin sce- 
nario, being an inverse polynomial, has no more then a 
finite number (the degree of the polynomial minus 2) of 
moments. At the same time, the corresponding asymp- 
totic pdf of the topological process, being of the Gibbs 
form oc exp(— $), may in principle admit all moments. 
If both mechanisms refer to a common asymptotic pdf, 
the latter being derivable in the Langevin approach, the 
previous statement is no longer valid, 

It turns out that the asymptotic behavior of a time- 
dependent pdf, in the semigroup (topological) modeling, 
may critically depend on the initial data choice, like e.g. 
the location of the initial pdf in the "potential landscape" 
of V* (x). The signature of such behavior is encoded in the 
functional form of a semigroup potential V(x), derived 
from the a priori chosen invariant pdf (x) , by means of 
above generalized reverse engineering procedure. If the 
effective (semigroup) potential obeys the requirements 
of [lCj . we may expect that a prescribed invariant den- 
sity is indeed approached in the large time asymptotic of 
the random process, irrespective of the initial (pdf) data 
choice. These requirements need to be verified for each 
specific guess about a functional form of the prescribed 
invariant pdf p*. 



we substitute y for fc) 



^n(y) = 



A,, 
±A, 



Ai(- 
, Ai( 



y), y > o 
-y), y<o, 



(Al) 



where n enumerates eigenvalues (and corresponding 
eigenf unctions). In Eq. (|A1[) . y n numbers the n-th zero 
of the function Ai(y) (for odd n), or of its derivative for 
even n (observe that zeros of both function Ai and its 
derivative lie on the negative semi- axis). A n is a normal- 
ization coefficient, determined by the standard identity 



A 



i }2 n {y)dy = !• 



(A2) 



Here we use simply ip 2 (rather then \ip\ 2 = ipip*) since 
the eigenfunctions are real. 

The method of construction of the above wave func- 
tions from initial Airy function is shown graphically in 
FigJS] We link either function (odd n - FigJ5]D) or its 
derivative (even n - Fig[9^,) in one of the zeros. In other 
words, the wave function of n-th state is created by shift- 
ing the Airy function on the positive half-axis to the right 
so that its n-th zero y n coincides with the origin (zero) 
of the coordinate system. After this step we continue a 
function to the negative half-axis either evenly (for even 
n) or oddly for odd n. The expression for A n is as follows 



A,- 



Ai 2 (-y n -y)dy + / Ai 2 (-y„ + y)dy 



(A3) 

The evaluation of integrals (|A3|) yields after some algebra 



1 

V2 



y n A?(-y n ) + Ai' 2 (-2/„) 



-1/2 



(A4) 

where prime means derivative with respect to an argu- 
ment. Further simplifications of Eq. (|A4[) are possible if 

we observe that for even n Ai 2 (— y n ) — and for odd n 
Ai 2 (-y„) = so that 



A n = 



[Ai'(-y„)V2] 1 , n is odd 
[Ai(-y n )s/2y^\ , n is even. 



(A5) 



APPENDIX A: SCHRODINGER EIGENVALUE 
PROBLEM FOR THE LINEAR POTENTIAL AND 
AIRY FUNCTION 



In the present Appendix we briefly recapitulate, ba- 
sically retrievable in the literature but not accessible in 
minute detail nor in a closed form, (see, e.g. [Hj]), a pro- 
cedure of construction of the eigenfunctions of equation 
(I27p . Symmetry arguments (see, e.g. [13) and an explicit 
solution of corresponding Schrodinger equation [25| lead 
to the following expression for the eigenfunctions (here 



Also, the corresponding (dimensional) energy eigenval- 
ues E n = Von can be found from a condition that in the 
p space all p n corresponding to zeros y n must be zero. In 
other words, x n = y n ( + a = or 



En 



\Vn\ 



2~ 



1/3 



(A6) 



Here we reflect the fact that zeros of the Airy function 
and its derivative are negative. 

Now we are in a position to write the explicit form of 
several first wave functions of Eq. (|2T)) . The ground state 
function is defined by the Eq. (|2"8|) . The first excited 
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FIG. 9: The idea of the wave functions construction; a- 
ground state (and even n), b - first excited state (and odd 
n) 



The third excited state (n — 3) reads 



Mv)=aJ A ^ m + y \ V> ° n (A9) 

w\yj ^ -Ai(-j/ 3 - y), y < o, v ; 

A, = Ui^-ya)^!' 1 , y 3 « 3.2482. 



The functions ([28]), (ATJ) - CM]) are plotted in Fig. [TO] 
The conformance with oscillation theorem (the n-th wave 
function of a discreet spectrum can have only n zeros on 
its domain) [25| is seen. 
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state (n = 1) has the form 



, Ai(-yi+y), y > 



Ar = [Ai'(-y x )v^ 
For the second excited state (n = 2) we get 
Ai(-y 2 + y), y > 



, yi « 2.3381. 



Mv) = a 2 , Ai( _ Ift ,_ y)iy<0) 



-4 2 



Ai(-y 2 )\/2^ 



y 2 « 3.2482. 



(A7) 



(A8) 



FIG. 10: Several wave functions for the linear modulus po- 
tential, Eq. (|27)> . 



APPENDIX B: THE APPROXIMATE 
EXPRESSIONS FOR ZEROES OF AIRY 
FUNCTION AND ITS DERIVATIVE 



For large negative y the function Ai(y) has the follow- 
ing asymptotic expansion 



Ai(x — > — oo) 



cos 



7T 

4 



2xy 



(-l) 3/4 (j) 1/4 5(-l)V*(£) 
V^r 48^ 



7/4 



■ sin 



7T 

4 



2x-v 



(Bl) 



Equating series (|B1I) to zero, we obtain the desired ana- cays at infinity much faster then that preceding the co- 
lytical expression for zeros of the Ai function. We observe sine. So, we simply equate to zero the argument of cosine 
that the coefficient before sin, proportional to x~ 7 / 4 , de- 
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function, getting 



2/3 



n + - 

4 



2/3 



(B2) 



It turns out that (|B2|) gives a fairly good approximation 
for zeros that begin at the lowest eigenvalue n — 1. We 
compare exact and approximate (|B2[) solutions: 



y§ xact 

1 y° xact 
n = 2 y° xact 



n 
n 



-2.3381, y appr = -2.32025, 
-4.0879, ?/ appr = -4.08181, 
-5.5206, y appr = -5.51716, 



n = 8 y^ xact = -11.0085, % ppl = -11.0077, 

.exact _ 1 1 noc ..appr ^ _H. 9 353.(B3) 



n = 9 y§ xact = -11.936, y aI 



The same asymptotic analysis can be performed for 
the derivative of Airy function. We end up with 



V 



(+) _ _ TL 



2/3 



1 

n + - 

4 



2/3 



(B4) 



A comparison of exact and approximate roots goes as 
follows 



n = y( +)cxact = -1.0188, y£ +)appr 
n=l y W° xact = -3.2482, y[ +)apw 
n = 2 y W cxact = -4.8201, j/< +)appr 



n = 8 y( +)cxact = -11.4751, J/' 



(+)appr 



n = 9 ?/g 



(+)cxact 



12.3848, y 9 



(+)appr 



-1.11546, 
-3.26163, 
-4.82632, 

-11.4762, 

-12.3857. 

(B5) 



(from momentum to coordinate space) Fourier transform 
is defined as f(x) = ± f(p)e~ lpx dp . 

For the ground state we proceed accordingly. We sub- 
stitute Eq. (|28p into the Fourier integral to obtain 



Zir 



(CI) 



Ai(-y -p)e ipx d P = / Ai(t)e- ix ( t+Vo ^dt, 



Ai(-y +p)e lpx dp 



M{t)e- lx{t+yo) dt, 



h+h= / Ai(t) 



o -ix(t+y ) _|_ gix(t+Vo) 



fif. 



Hence, the Fourier image of the ground state wave 
function is determined by the above equation (|29| . For 
the higher even n the above method yields 



'0even(*^) — 

7T 



For odd states 



Ai(i) cos ae(* + y n )dt, (C2) 



^odd(p) — 7T~"(^2 — -Tl)j 
Z7T 



/ 2 



Ai(-2M-p)e^: 
Ai(-y„ + p)e 4p:c dp = 



(C3) 

DO 

ki{t)e- ipi - t+v - ] dt, 

Vn 
JO 

Ai(t)e ix( - t+Vn Ut, 



h-h = A n I Ai(t) 



e ip(t+y n ) _ e -ip(t+y„ 



dt. 



In other words, the Fourier images of odd wave functions 
have the form 



APPENDIX C: FOURIER IMAGES OF THE 
WAVE FUNCTIONS 



V'odd(p) = *4r / Ai W sin P( f + yo)dt. (C4) 

y n 



Since our departure point in Section III D was 
the Schrodinger-type eigenvalue problem in momentum 
space, it turns out to be useful to discuss Fourier images 
of the above eigenfunctions. We recall that the inverse 



Here, A n are determined by Eqs. (|A5[) . It is seen that 
odd Fourier images are imaginary odd functions. This 
(imaginary coefficient) does not affect the physical mean- 
ing of IV'oddCp)! 2 which is a probability density. 
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